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Abstract 

We discuss the implementation of arbitrary precision composite pulses developed 
using the methods of Brown et al. [Phys. Rev. A 70 (2004) 052318]. We give explicit 
results for pulse sequences designed to tackle both the simple case of pulse length 
errors and for the more complex case of off-resonance errors. The results are de- 
veloped in the context of NMR quantum computation, but could be applied more 
widely. 
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1 Introduction 



Composite pulses pQ have long played an important role in many NMR exper- 
iments, allowing the effects of systematic errors to be reduced. More recently 
they have been applied in quantum computing [2], including both magnetic 
resonance experiments and other implementations [3f4l5ll6] . Composite pulses 
developed for quantum computing differ from more conventional NMR ap- 
proaches in two key ways. Firstly they must perform the desired rotation 
whatever the starting state of the system, so that they act as general rotors; 
in NMR such pulses are sometimes called type A composite pulses PQ, and 
have the advantage that they can be inserted into any part of a pulse sequence 
without the need for careful analysis. Secondly, they are usually designed to 
give extremely accurate rotations in the presence of small errors, rather than 
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to give reasonable rotations in the presence of large errors. Despite these dif- 
ferences, however, composite pulses designed for quantum computation can 
be used in conventional NMR experiments, and studying them can give some 
insight into more conventional approaches. 

Type A composite pulses are most commonly used to tackle pulse length 
errors, which occur when the strength of the driving field differs from its 
nominal value, for example as a result of inhomogeneity, and off-resonance 
errors, which occur when the frequency of the driving field is not quite in 
resonance with the transition of interest. For pulse length errors the BB1 
sequence j7], originally developed by Wimperis, has proved highly successful. 
Designing good type A pulses for off-resonance errors has proved more difficult, 
with the CORPSE sequence [HE], based on an early numeric result by Tycko 
[8], being perhaps the most satisfactory. 

More recently Brown et al. have described a general method PfTU] for gener- 
ating arbitrary precision composite pulses, which they claim can be applied to 
tackle errors of arbitrary kinds. They do not, however, give explicit solutions 
for more than a small number of cases. We will show that it is simple to apply 
their methods to generate pulses resistant to pulse length errors, but it is more 
difficult to tackle off-resonance effects. 



2 Pulse length errors 

We consider rotations about axes in the xy plane, for which an ideal rotation 
is described by the propagator 

U (9, (f>) = exp[— i 9{a x cos0 + a y sin 0)/2] (1) 

where 9, the rotation angle, depends on the strength of the resonant driving 
field and the time for which it is applied, and 0, the phase, depends on the 
phase of the field with respect to some appropriate reference. Following Brown 
et al. [9] we describe our rotations in terms of the Pauli matrices, which are 
trivially related to the corresponding single spin product operator terms 
The description above is really only appropriate when 9 > 0, but also works 
for negative angles as 

U(-9,^ = U(9^ + n) = U-\9, ( f ) ). (2) 

In the presence of errors, however, this extension may not be appropriate and 
it is necessary to proceed with some caution. 
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Pulse length errors occur when the strength of the field is higher or lower than 
its nominal value, so that all rotation angles are systematically wrong by some 
constant fraction e. It is convenient to write 

v{e,<i>) = u{e[i + eU) (3) 



in such cases, and a series expansion in the error e gives 

V{9,<f>) = U(e,<f>) + 0(e) (4) 



so that the propagator has first order errors in e. Note that in this paper 
we write all our composite pulses as sequences of propagators, so that the 
order of pulses runs from right to left. We classify pulse sequences according 
to n, the order of error in the propagator, but will occasionally refer to the 
corresponding propagator fidelity 

f= 1 -Tr(VU- 1 ) = l-0(e') (5) 



in which errors appear to order 2n. For pulse length errors 

V(M + 7r) = ^ -1 (M) (6) 



which will prove very useful throughout this section. 

Several composite pulse methods for correcting such pulse length errors exist, 
most notably BB1 [7j. This replaces a single pulse with a sequence of pulses 
such that 

V(tt, <P a )V(27f, <f) b )V(7T, <f> a )V(e, 0) = U{0, 0) + 0(e 3 ) (7) 



where 0& = 30 a and a = arccos(— 0/Att), to give a propagator with third order 
errors in e. This pulse sequence was discovered using geometric arguments, and 
turns out to be remarkably effective in practice [4"][T2] . 

Brown et al. replaced previous methods of finding composite pulses, based 
on intuitions or special forms, by a systematic procedure j9]. We begin by 
describing in detail how this procedure can be used to generate a series of 
composite pulses to correct pulse length errors. As before we only consider 
target rotations with phase angles of zero, as more general rotations in the 
xy plane are trivially derivable from these by offsetting the phase of all pulses 
appropriately. These pulses can also be used to design sequences for robust 
evolution under J-couplings [T3] . 
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2.1 Isolating the error 



The first step is to isolate the error part of a pulse sequence by calculating 

7/9 

A x = V(9, 0) U-\e, 0) = 1 - I e a x + 0(e 2 ) ({ 



where A\ indicates a first-order error in a sequence correct to zero-order. Since 
A- 1 V = UV- 1 V = U (9) 

if we can generate A^ 1 exactly then we can convert an erroneous rotation into 
a correct one. More realistically, if we can generate A x l correct to first order 
then we can use this to cancel out the first order error in the original pulse. 
To do this it is useful to note that 

A x = U(e9, 0) + 0(e 2 ) (10) 
so that A^ 1 is approximately equal to an x-rotation with angle — 9e. 

2.2 Generating a pure error term 



To correct the first-order error it is necessary to generate a matching rotation 
whose angle is directly proportional to the error fraction e. A pure error term 
of this kind is most easily generated by noting that 

K(2tt,0) = -l + ineo- x + 0(e 2 ) (11) 



has the correct general form (the global phase of —1 can be ignored as usual). 
However while the error has the right form it has the wrong magnitude. This 
can be fixed by using two rotations with different phases [9], as 

Xx(0) = V(2n,(p)V(2TT, -0) = 1 -i2ix cos(<p)e a x + 0(e 2 ) (12) 



and so the size of the error can be scaled. In particular, solving 

-2ttcos(0i) =9/2 (13) 

to get 

0i = ±arccos(-6»/47r) (14) 
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allows the first order error to be corrected (note that the sign of the error term 
must be reversed as we are seeking to approximate A^ 1 rather than A\). This 
gives 

X 1 ((ft 1 )V(6,0) = U(6,0) + O(e 2 ) (15) 

as a pulse sequence which is correct to second order. Interestingly the key 
phase angle 4>i has the same size as (ft a in the BB1 sequence. The sign of 4>i 
(and (ft a ) seems arbitrary, but for definiteness we choose the positive value. 
This point will be examined in more detail later. 

The size and scaling of the errors permits a value of (ft± to be found for any 
choice of 6, but if larger error terms are needed they can be achieved by simply 
repeating the sequence of 2ir rotations, thus doubling the error. In passing we 
note that although the approximate form given in Eq. [T2] is independent of 
the sign of (ft, and thus of the order of the two 2tt pulses, the exact form of 
X\ does in fact depend on the sign of (ft, and it is therefore necessary to use a 
consistent order. This will become important below. 

2.3 Treating the second order error 

The process described above can then be repeated to isolate the second order 
error 

A 2 = X 1 ( ( ft 1 )V(6,0)U- 1 (6,0) (16) 
with the result 

A 2 = l e 2 a z + 0(e 3 ). (17) 

o 

To correct this will require an error term dependent on e 2 , rather than on e, 
and directed along the z-axis rather than the x-axis. This can be achieved 
by using the properties of group commutators. As these properties play a key 
role in the proof of the Solovay-Kitaev theorem [TJ] Brown et al. refer to the 
results as SK pulse sequences, but beyond the role of group commutators there 
is no need to undertand the Solovay-Kitaev theorem to see how their methods 
work. The key result 

exp(-iAe l ) exp(-iBe m ) exp(iAe') exp(ifle m ) 

= exp([A B]e l+m ) + 0(e l+m+1 ) (18) 
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shows that two different first order pure error terms can be combined to make 
a single second order pure error term as long as their directions are chosen 
properly and inverses are available for each error term. This restriction will 
not be important for pulse length errors but will be more problematic for 
off-resonance errors. 

As the commutator [a x , a y ] = —2ia z the desired second order error term along 
z can be generated from x and y terms. The x term can be generated as 
before, Eq. [TjJ and the equivalent y term Y\ can be generated in the same way 
by shifting the phase of both pulses by ir/2. (Brown et al. in fact described 
an alternative method [9 J for implementing Y±, but this alternative is more 
complex.) Inverse terms can be generated by reversing the order of the two 2tt 
pulses and shifting their phases by n. Note that it is necessary to reverse the 
sequence of pulses even though Eq. [T2] appears to be independent of this, as 
the two alternatives differ to second order in e; it is not sufficient to generate 
an inverse which only accurate to first order. Hence 

ZM = X 1 (<l>)Y l ( ( j ) )X^( ( j ) )Y 1 - 1 ( ( f ) ) = 1 - z8tt 2 cos 2 (0) e 2 a z + 0(e 3 ) (19) 



allows a second order z error of some desired size to be generated. If an error of 
the opposite sign is needed then the X\ and Y\ sequences can be interchanged 
to give 

Z' 2 {4>) = Y 1 ((f>)X 1 ( < f>)Y 1 -\(f>)X^ 1 ((l>) = 1 + *87r 2 cos 2 (0)eV z + 0(e 3 ). (20) 



As before we can solve 

8tt 2 cos 2 (0 2 ) = eVl6n 2 -9 2 /8 (21) 

to remove the second order error, giving 

(22) 



S>2 = ± arccos 



•57T 



where the signs may be chosen independently, and we choose initially to take 
both positive signs. Thus 

z' 2 ( ( p 2 )x 1 ( ( f )1 )v(e, o) = u(e, o) + o(e 3 ) (23) 

is our desired pulse sequence, correct to third order. The explicit expansion 
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V(2ir, tt/2 + 2 )y(27r, tt/2 - (p 2 )V(2n, <p 2 )V(2n, -02) 
^(2vr, 3vr/2 - 2 )V(2vr, 3vr/2 + 2 )V(2vr, vr - 02)\/(2vr, vr + 2 ) 

l/(27r,0O^(27T,-01)V(0,O) (24) 



shows that this contains a total of ten correction pulses in addition to the 
main pulse. 

2.4 Rotating and redividing the error 

The sequence described above is not in fact the sequence originally described 
by Brown et al. Their approach [9] is instead based on generating an X 2 error 
correction term using appropriate Y\ and Z\ sequences and then rotating this 
onto the z axis. In the absence of systematic errors, such rotations are easily 
performed. For example the identity 

[/(vr/2, 3tt/2)£7(0, 0)U(n/2, tt/2) = exp(-z0 a z /2) (25) 

allows a z rotation to be generated from x and y rotations, a composite Z-pulse 
[j~6] . In the presence of systematic errors it is necessary to proceed with more 
caution, as the errors in different pulses will combine in a complex manner. 
However, in the presence of pulse length errors it is possible to use imperfect 
pulses to rotate pure error terms, as 

V(tt/2, 3ir/2)X n V(n/2, tt/2) = Z n + 0(e n+1 ) (26) 

for any error order n, and similarly for other error terms. This approach re- 
quires the ability to generate accurate inverse operators, and while this is easy 
for pulse length errors it can be tricky in other cases. 

For pulse-length errors we can implement the second order correction sequence 
using 

Z 2 (0 2 ) w V(tt/2, 3n/2)X 2 (4> 2 )V(n/2, tt/2) (27) 

X^(0 2 ) = Z 1 (0OF a (0i)Zr 1 (0Oyr 1 (0O (28) 
Zi{fa) « V{tt/2, 3vr/2)X 1 (0 2 )V(7r/2, tt/2). (29) 

This alternative sequence gives identical performance at second order (com- 
plete correction of errors) but differs in its third order behaviour, as we shall 
see later. 

We can also consider many other possibilities: firstly we can instead generate 
Z' 2 from any of Y 2 \ X 2 , or X 2 , secondly we can use alternative rotations to 
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generate Z\\ thirdly we can use the negative sign for 0i in the first order cor- 
rection sequence (in which case the second order error changes sign). Beyond 
these possibilities, built around rotating the error, we can also choose how 
to divide up the relative contributions to the second order error term arising 
from the two first order terms. For example we can write 



Z 2 (a,(3)=X 1 (a)Y 1 ((3)X^\a)Y 1 -\P) 

= 1 - i 8tt 2 cos(a) cos(/3) e 2 o z + 0(e 3 ) (30) 

and control the size of the second order error term by varying a and f3. More 
simply still, we can use the fact that V(2tt, 0) gives an unsealed pure error 
term along x to use the form 

V(2tt 1 0)Y 1 (/3)V(2tv 1 7i)Y 1 - 1 ((3) = 1 - iAn 2 cos/?e 2 a z + 0(e 3 ) (31) 



which only has six pulses, rather than the usual eight, and obtain the correct 
error term by choosing 



(3 = ± arccos 



32^ 



(32) 



2.5 Choosing between sequences 



Given this plethora of subtly different sequences it is reasonable to ask which 
is the best. In some sense all second-order pulse sequences are equally good, as 
they all suppress errors to the same order, but it is possible to choose between 
them either by considering higher order errors, or by considering sensitivity 
to other types of error [I]. 

Here we adopt the first approach, choosing to minimize the size of the third 
order error term, and initially specializing to the case of 180° pulses, so that 
9 = 7i. The smallest error term we have so far been able to find occurs when 
using equations [27] to [29] to create the second order term, and taking positive 
signs throughout equations HH and [22] There is no obvious reason why this 
choice is best, but it does give a third order error more than 15 times smaller 
than some other alternatives. 

Interestingly, the second best choice we have found is the BB1 sequence, which 
for the case 9 = n has an error only about 10% larger than the best sequence. 
Furthermore for BB1 the size of the third order error term scales approximately 
linearly with 9, while the behavior of the "best" sequence, described above, 
is more complex. Thus for most flip angles (specifically, 9 < 168°) BB1 is the 
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best second order sequence known. For the case of 180° pulses it has a fidelity 



2.6 Third order errors 

We can of course correct the third order error term in much the same way as 
the first and second order errors. As pointed out by Brown et al. there is no 
need to use a fully systematic approach of correcting error orders in sequence; 
instead we can begin with BB1 and correct the third order error term. 

The third order error term for BB1 (and indeed for all the other sequences 
considered above) lies in the xy plane, with the size and position depending 
on the value of 9. Here we do not give complete results, but simply sketch 
a partial solution. There are many possibilities for generating a third order 
error, but one simple example is 



XM = Y 1 ( ( j ) )Z 2 ( ( f ) )Y 1 -\4 > )Z^^) = 1 - z32tt 3 cos 3 (0) e 3 a x + 0(e 4 ) (34) 



where Z 2 ((f>) is generated from Xi(<p) and Yi(<f>) as before, equation [191 an d 
Z 2 l = Z' 2 , so that this sequence only requires x and y rotations. The error 
term can then be rotated into the correct position, either by composite z 
rotations or, more simply, by shifting the phases of all the pulses in the X 3 
term. 

As before we need to choose a value of 03 to cancel the third order error, but 
as the calculations become very complicated we here consider only the special 
case of 180° pulses, 9 — tt. Even in this case the analytic result is complicated, 
and so we simply give the numerical value, 03 ~ 73.1°; the required phase 
shift can also be calculated as approximately —1.6°. 



3 Off-resonance errors 

The treatment of off-resonance errors is superficially similar but much more 
difficult in practice. The fundamental problem is that, unlike the case of pulse 
length errors, it is not possible to generate perfect inverses of arbitrary rota- 
tions in the presence of off-resonance errors. 

Off-resonance errors occur when the frequency of the driving field is not quite 
in resonance with the transition of interest, so that all rotations occur around 
a tilted axis. They can be parameterized in terms of the off-resonance fraction 




(33) 
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/, equal to the ratio of the frequency error and the driving frequency. The 
basic rotation in then 

V(6, (f)) = exp[-i 9(a x cos + a y sin + / a z )/2] = U(0, (f)) + 0(f) (35) 



which has first order errors in /. Unlike the case of pulse- length errors this 
definition should only be used for positive values of 9. 

The CORPSE family of sequences [3|4"] for correcting off-resonance errors uses 
the three pulse sequence 

c(6, 0) = v(e c , <p)v(e b , + tt) v(e a , <p) = u(e, <p) + o(f) (36) 



where 



9 a = n a 2n + 9/2 - arcsin[sin(#/2)/2] (37) 
6 = n 6 27r-2axcsin[sin(0/2)/2] (38) 
9 c = n c 2-n + 9/2 - arcsin[sin(#/2)/2] (39) 

and the size of the second order error term depends on the values chosen 
for the integers n a , and n c . The smallest errors are seen for the original 
CORPSE sequence, which has n a = = 1 and n c = 0; for the case of 180° 
pulses the fidelity is 

F=l-*±±M,< + 0( r). (40) 



Short-CORPSE, defined by n a = n c = and n& = 1 is the shortest possible 
sequence but has a much larger error term pQ. 

The CORPSE family will play a key role in the following sections, not simply 
because it provides a pulse sequence with no first order errors, but mostly 
because it provides a route to sufficiently accurate inverse propagators. In 
the presence of off-resonance errors V(9, <p + tc) is not an accurate inverse for 
V(9,<f>) as 

V(9,7r)V(9,0) = l + O(f). (41) 

The corresponding CORPSE pulses perform much better, 

C(9,n)C(9,0) = l + O(f), (42) 
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and provide sufficiently accurate inverses to allow error terms to be rotated as 
for pulse length errors. As before the size of the third order error term depends 
on the exact choice of sequence, but is now smallest for short-CORPSE. 



3. 1 Correcting first order errors 

We now explore the systematic correction of error orders using the methods 
previously described. The first order error can be isolated as before, 

A x = V (9, Op' 1 (9,0) 

= l-i sin(0)/2 fa z + i sin 2 (#/2) f a y + 0(f) (43) 

and in general lies in the yz plane. Tunable pure error terms can be created 
either by using the form given by Brown et al. [9], 

B x (4>) = V(<f>, 0)V(2cf>, ir)V(4>, 0) = 1 - i 2 sin(0) / a z + 0(f), (44) 

or the alternative form 

Y[(<P) = V(n, <t>)V(ir, vr + <f>)V(n, -</>)V(ir, n - <f>) 

= l + iA C os(<P)fo y + 0(f). (45) 



Designing a sequence for the case 9 = it is easy as the error lies solely along 
the y axis in this case, and so 

Yl(^)V(n,0) = U(n,0) + O(f) (46) 



withe the choice <\>\ = arccos(— 1/4) rs 104.5°. Interestingly, the key phase 
angle in this sequence turns out to be the same as that used in a BB1 pulse 
with = 7i. The size of the second order error term is significantly larger than 
for CORPSE (and somewhat larger than for short-CORPSE), with a fidelity 

?=l-^ff + 0(f). (47) 



but this sequence does have the relative simplicity of being constructed solely 
from 180° rotations, albeit with complicated phases. 

Designing a sequence for other values of 9 is, however, much trickier. In the 
general case the error does not lie along a principal axis, and so it might seem 
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that we should rotate one of the two pure error terms. This cannot be done 
directly using simple rotations, as accurate inverse propagators are required 
to rotate error terms. It could be achieved using CORPSE pulses, but this is 
not sensible as CORPSE is already correct to first order. 

An alternative approach is to note that pure error sequences can simply be 
combined, and so build up a tilted error by combining the z and y error 
sequences, 

B 1 {<t>T)Yl{4>l) = 1 - i 2 sin(0f ) / a z + 1 4 cos(0?) f a y + 0(f) (48) 

with any cross terms between the two parts if the pulse sequence being swal- 
lowed up into the 0(f 2 ) term. Choosing 0f = arccosf— sin 2 (6>/2)/4] and 
<fil = — arcsin[sin(#)/4] allows first order off- resonance errors to be suppressed 
in the general case. However, the size of the second order error term remains 
significantly larger than for CORPSE, and this general sequence does not have 
the simplicity seen for the special case of 180° pulses. Thus CORPSE remains 
the best currently known type A sequence correct to first order in the presence 
of off-resonance errors. 



3.2 Correcting higher order errors 

The systematic approach can, more sensibly, be used to correct higher order 
errors. As the equations for arbitrary values of 6 become extremely compli- 
cated we will again limit ourselves to the special case 9 — n. We extend our 
definitions 

X x (0) = V(n, tt/2 - 4>)V(tt, 3tt/2 - 4>)V(tt, tt/2 + (/>)V(tt, 3tt/2 + 0) 

= l-i4cos(</>)f<T x + 0(f 2 ) (49) 

1^(0) = V(ir, ir - <j>)V(ir, -<f>)V(ir, tv + <j>)V(ir, 0) 

= l-*4cos(0)/^ + O(/ 2 ) (50) 

X[{4>) = V(ir, 3tt/2 + 4>)V(ir, tt/2 + (p)V(n, 3tt/2 - </>)V(ir, tt/2 - 0) 

= 1 + 24cos(0)/ ( t ;c + O(/ 2 ) (51) 

(52) 

and note that X[((f>) ~ Xf 1 (0). The approximation is good enough that these 
terms can be used to build higher order propagators; in particular 

Z 2 (0) = X 1 (0)y 1 (0)X((0)y i '(0) = 1 -^32cos 2 (0)/ 2 a 2 (53) 
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and 

Z' 2 {<j>) = Y^X^Yl^X'M = 1 + *32cos 2 (0) fa z . 



(54) 



The second order error can be isolated as usual 



A 2 = Y(( ( j> 1 )V(7r,0)U- 1 (n,0) 

= l-i VlE/2 f<r z -i tt/4 f a x + 0(f) 



(55) 



and lies in the xz plane with magnitude V60 + 7r 2 /4. 



One approach to correcting this is by using CORPSE pulses to rotate an 
appropriate Z 2 error around the y axis to get the final sequence 



C(V> 2 , 7T/2)Z^ 2 )C(^ 2 , 37r/2)r i / ( ( /. 1 )\/(7r, 0) = U(n, 0) + 0(f) (56) 



This approach can, of course, be generalized to other angles and higher orders, 
but the resulting algebra is very complex. For simplicity it is possible to check 
results using ideal rotations in place of CORPSE based rotations: this will 
give the right result for error terms which are completely suppressed, but the 
wrong values for higher order errors. 

Alternatively, we can construct the tilted error term out of a combination of 
z and x errors as in the previous section. We begin by noting that Bi(3ir/2) 
provides a good inverse for the pure error term Bi(n/2), and that this allows 
us to construct a second order x error using 



with 




(57) 



and 




(58) 



X 2 (<j>) = Y, (0)B 1 (7i/2)Y((4>)B 1 (3tt/2) = 1-U6cos(0) fa, 



(59) 



This can be combined with a z error to get the pulse sequence 



X^fflZttfflYlMVfaO) = U(n,0) + O(f) 



(60) 
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where 0f = arccos(-7r/64) « 92.8° and 0* = arccos(v / 15/8) w 75.8°. This 
sequence has the advantage over the CORPSE based approach of requiring 
only 90° and 180° pulses. 

3.3 Time symmetric sequences 

In passing we consider the use of time-symmetry in composite pulse sequences 
[I]. In the presence of off-resonance errors time symmetric composite pulses 
have fidelities which are even functions of the off- resonance fraction / |15j . 
and so give the same fidelity for +/ and — /, although the details of the error 
may differ. Although such symmetric fidelities have no advantage in principle, 
the results are certainly easier to interpret. 

As an example we consider a time symmetric version of the pulse sequence 
to correct first order errors arising from off-resonance effects in a 180° pulse, 
which takes the form 



where <p' x = arccos(— 1/8) ~ 97.2°. The fidelity of this sequence has no fifth or- 
der term, unlike that of the previous version, equation H6| but as both fidelities 
are dominated by fourth order errors this is largely a cosmetic improvement. 

3.4 Simultaneous errors 

So far we have only considered the effects of pulse length errors and off- 
resonance errors in isolation, while in real physical systems both sorts of error 
can occur simultaneously. It is generally difficult to find pulse sequences which 
suppress both sorts of error simultaneously, but it is still important to consider 
whether insensitivity to one type of error is obtained at the cost of increased 
sensitivity to other types of errors [I]. 

As noted previously [I], the response of the time symmetric version of BB1 
to off-resonance errors is very similar to that of a simple pulse. This occurs 
because in the absence of pulse length errors the correction sequence 



has no first order terms arising from off-resonance errors, and so does not 
contribute significantly to the total error. In the same way in the absence of 




(61) 



V{ir, <Pi)V(2n, 30 1 )F(7r, 00 = 1 + O(f) 



(62) 



14 



off-resonance errors the correction sequence 

V(7T, 4> X )V{k, 7T + 1 )y(7T, -0 1 )y(7T, 7T - = 1 (63) 

has no error terms arising from pulse length errors at all. Thus these two 
sequences can be combined, giving the composite 180° pulse 

V(tt, faWn, vr + 1 )y( 7 r, -<h)V(ir, vr - fa) 
V{ir,fa)V{2<K,?>fa)V{<K,fa)V{<K^) 

with 0! = arccos(— 1/4). This has a fidelity 

15 4 5tt 6 6 169tt 2 22 

y- = 1 f e / e + higher terms 

8 1024 32 J y 

so in the absence of off-resonance errors the correction of pulse length errors is 
identical to a BB1 sequence, and in the absence of pulse length errors the cor- 
rection of off-resonance errors is even better than the simple pulse, equation l4~6l 
This pulse can correct well for either pulse length errors or off-resonance er- 
rors; in the presence of simultaneous errors the performance is not so good, 
but is still much better than a simple 180° pulse. 



4 Summary 

The methods of Brown et al. can indeed be used to derive arbitrary precision 
composite pulses, but the process can be somewhat complicated. For the case 
of pulse length errors the situation is simple, as it is possible both to gener- 
ate a wide range of pure error terms and their inverses, and to rotate these 
terms using uncorrected pulses. The case of off-resonance errors is much more 
complicated because the difficulty of generating accurate inverses of incorrect 
rotations means that the most direct approach cannot be used. 
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